Photon-assisted electron transport in graphene 
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Photon-assisted electron transport in ballistic graphene is analyzed using scattering theory. We 
show that the presence of an ac signal (applied to a gate electrode in a region of the system) 
has interesting consequences on electron transport in graphene, where the low energy dynamics is 
described by the Dirac equation. In particular, such a setup describes a feasible way to probe energy 
dependent transmission in graphene. This is of substantial interest because the energy dependence 
of transmission in mesoscopic graphene is the basis of many peculiar transport phenomena proposed 
in the recent literature. Furthermore, we discuss the relevance of our analysis of ac transport in 
graphene to the observability of zitterbewegung of electrons that behave as relativistic particles (but 
with a lower effective speed of light). 

PACS numbers: 73.23.Ad,73.23.-b,73.63.-b,03.65.Pm 



I. INTRODUCTION 

Since the discovery of an anomalous quantum Hall 
effect in graphene [l|, |2|, promising possibilities to ob- 
serve quantum dynamics of Dirac fermions in such sys- 
tems have been proposed. The most prominent one 
is the quantum-limited conductivity (of order e 2 /h), 
where measurements p], Q and theoretical predictions 
[E [J IE IE El are still inconsistent with each other. 
Other examples of unusual quantum transport phenom- 
ena in mesoscopic graphene are a maximum Fano fac- 
tor of | [9|, selective transmission of Dirac electrons 
throug h n-p junctions |ld| . the phenomenon of Klein tun- 
nelling [ri|. and transport phenomena at interfaces be- 
tween graphene and a superconductor 0, 0, Il5| . 
Recently, the effect on the longitudinal and Hall conduc- 
tivity of an applied microwave signal has been analyzed 
in bulk graphene |16| . In many of these works, interesting 
phenomena arise because transport through graphene is 
in general energy dependent. Therefore, it is desirable to 
directly probe the energy dependence of transport. We 
show that the Tien-Gordon problem |l7L Il8| of photon- 
assisted electron transport for Dirac electrons is a power- 
ful tool to quantitatively probe energy dependent trans- 
mission in graphene. The idea is to apply an ac signal to 
a region of graphene and a bias with respect to a neigh- 
boring region to allow electron transport from the one to 
the other (see Fig.HJ 

In the dc limit, electron transport then allows to di- 
rectly determine the energy dependent transmission co- 
efficients of the underlying scattering problem. Addi- 
tionally, we find that resonance phenomena (sharp steps 
in dG/dV, where G = dl/dV is the differential conduc- 
tance) arise if the applied bias V equals multiples of the 
applied ac frequency uj. These resonance phenomena are 
due to the vanishing of propagating modes directly at 
the Dirac point (the point in the spectrum of graphene, 
where the valence band and the conduction band touch 
each other). This implies an interesting application of 



graphene: It can be used as a spectrometer for high fre- 
quency noise - similar to the superconductor-insulator- 
superconductor (SIS) junction in Ref. [2(J with the ad- 
vantage that there is no frequency limit. 

Another motivation to look at the Tien-Gordon prob- 
lem in graphene is its relevance to the observability of an 
interesting and unobserved phenomenon, present for free 
Dirac fermions but absent for free Schrodinger fermions, 
called zitterbewegung (ZB) |2l|,|22j. ZB manifests itself, 
for instance, in a time-dependent oscillation of the posi- 
tion operator of a Dirac electron in the Heiscnberg pic- 
ture. A pioneering attempt to describe an experimental 
way to observe ZB in III-V semiconductor quantum wells 
has been proposed in Ref. 1231 . Others have adopted this 
idea to carbon nanotubes ]24|, spintronic, graphene, and 
superconducting systems |25j. Another situation where 
non-relativistic electrons experience ZB is the presence 
of an external periodic potential (2(| • In fact, there 
is a wide class of Hamiltonians (where the correspond- 
ing wave function is a spinor) which all exhibit ZB |27| . 
Nevertheless - although almost omnipresent - the phe- 
nomenon has never been observed in nature. Therefore, 
it is desirable to propose an unambiguous signature of 
ZB in a particular setup. As we show below, the condi- 
tions to observe ZB can be easily created using the Tien- 
Gordon effect, which results in an oscillating current due 
to the existence of a wavefunction that is a superposition 
of different energies. However, an unambiguous detection 
of ZB in this way is complicated, since also for the free 
Schrodinger equation an oscillating current is generated 
when the wavefunction is a superposition of different en- 
ergies. Additional studies are needed to single out effects 
of ZB and ensure their unambiguous observation. 

The paper is organized as follows: In Sec. [HI we in- 
troduce the setup under consideration and the model to 
describe it. The photon- assisted current is analyzed in 
Sec. IIIII Afterwards, in Sec. IIVI we discuss the rele- 
vance of the previous analysis to the observability of ZB 
in graphene. Finally, we conclude in Sec. 
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FIG. 1: Schematic of the setup. Two regions of graphene 
(called Gr(ac) and Gr (in)) are voltage-biased with respect to 
each other with a bias V. A gate electrode (grey-shaded area) 
shifts the Dirac point of region Gr(ac) away from zero energy. 
Additionally a small ac signal is applied to the gate to gen- 
erated side band transitions. The applied gate voltage V sa ,(t) 
has to be chosen in such a way that the (time-dependent) 
chemical potential in region Gr(ac) is given by Eq. 0. 



II. SETUP AND MODEL 

The setup under consideration is illustrated in Fig. ^ 
It contains two regions of graphene called Gr(ac) and 
Gr (in). A gate electrode in region Gr (ac) shifts the Dirac 
points of the two regions with respect to each other by 
an amount eVg. Additionally, we apply a small ac signal 
to the gate electrode. 

An ideal sheet of graphene (in the absence of K—K' 
inter-valley mixing) can be described in the effective 
mass approximation by a two-dimensional Dirac equa- 
tion (DE) for a two-component wave-function envelope 
\& = ^2) (subscript 1, 2 refers to pseudospins, whose 
origin can be traced to the presence of two carbon sub- 
lattices) 
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where v is the Fermi velocity. The dimensions of the sam- 
ple are assumed to be large enough such that boundary 
effects can be neglected. It is, however, straightforward 
to include boundary effects along the lines of Ref. 0. 
The index e/h refers to electron- like (energy e > with 
respect to the Dirac point) and hole-like (energy e < 
with respect to the Dirac point) solutions to the DE. The 
chemical potential of the two regions Gr(ac) (x < 0) and 
Gr(in) (x > 0) is 



eV s 



eV ac cos(wi) , 
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Gr (ac) and transmitted to Gr(in). The solution to the 
wave equation for x < (in region Gr (ac)) may then 
be written as 
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(We have dropped the argument x of the wave function 
in the latter equation.) In Eq. Q, J m is the mth order 
Bessel function. ^,From now on, we focus without loss 
of generality on electron-like solutions only and drop the 
index e. This is justified because, in ballistic transport 
from one region of graphene to another region, particles 
with a fixed energy s (which can be either positive or 
negative with respect to the Dirac point) are transmitted 
and their energy is conserved in the absence of inelastic 
scattering. 

The electron-like plane wave solutions of Eq. Q can 
be written as linear combinations of the basis states 
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where 



a ac = arcsm 



hvq 
e + eVsJ ' 



(5) 
(6) 

(7) 



q is the transversal momentum, and fc ac = (e + 
eVs) cosa ac /(hv) the longitudinal momentum in region 
Gr (ac). Likewise, electron-like plane wave solutions of 
Eq. Q in region Gr (in) can be written as linear combi- 
nations of the basis states 
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with a>i n = &rcsui(hvq/e) and k- ln = (e/Kv) cos a- ln . Now, 
we solve the transmission problem from region Gr (ac) to 
region Gr (in). An incoming wave function from region 
Gr (ac) is given by 



m— — 00 



(2) The reflected wave function in region Gr(ac) reads 



The potential fi(x) is chosen such that a solution to the ^^ ac \x,t) 
DE with components of different energies is generated in 
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where r m is the energy-dependent reflection coefficient. 
Furthermore, a transmitted wave function in region 
Gr (in) can be written as 



m— — oo 
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where 



with 
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hvq \ 
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(10) 
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(12) 



and fcin,m = (e + Timuj) cos a m , m / (fiv) . In Eq. I|I0(I . t m is 
the energy-dependent transmission coefficient. In order 
to determine t m and r m , we need to match wave functions 
at x — 0, namely 

*< ac) (x = 0, y , t) + *(. ac ) (x = 0, y, t) = *<i n) (x = 0,y,t). 
The solutions to the resulting set of equations are 



1 _|_ gi(a ac +Qin, m ) ' 
-i(Q ac -Q in , m )/2 /COSain.n: 



(13) 



1 + 
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1 _|_ gi(aac+ain, m ) ' 



where a ac is given by Eq. Q) and a; njm by Eq. (|12l) . It 
is easy to verify that unitarity holds \r m \ 2 + \t m \ 2 = 1. 
Furthermore, Eq. I|13fl shows that if the angle of incidence 
is zero (q = 0) then all transmission coefficients are 1 and 
all reflection coefficients vanish. This is known as Klein 
tunnelling in relativistic quantum dynamics [TT| . 



III. PHOTON-ASSISTED CURRENT 

In order to determine the transmitted current, we need 
to calculate the current density operator in n-direction 



J x (x,t) = ev%?*(x,t)<7 x V{x,t), 



(14) 



integrate over a cross section in y-direction, and over an- 
gles of incidence. In the following, we assume that a dc 
bias V is applied between regions Gr(ac) and Gr(in) and 
that ksT is the lowest of all energy scales. The average 
current is then given by 
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where W is the width of the sample and g max 
Ml II {\eV + hmu\/hv, \eV + eV s \/hv\ is the upper bound 
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FIG. 2: (Color online) dG/dV is plotted in units of 
8e 3 W/(vh 2 ) for a fixed value of eVs/hu = 100 and differ- 
ent values of eV ac /Tiu>. Sharp steps appear when the applied 
bias equals multiples of the ac frequency. The size of the steps 
is non-universal. It depends on the weight of the different side 
bands given by the Bessel functions in Eq. 1161 . For the given 
values of eV ac /huj, it is sufficient to include side bands up to 
\m\ — 4 to calculate dG/dV to a very high accuracy. 



of the transversal momentum of propagating modes. The 
bias is applied such that the Fermi function (at zero tem- 
perature) in region Gr (ac) reads f ac (E) — 6(eV — E) 
and, in region Gr(in), fi n (E) = 6{—E), where 8{x) is the 
Heaviside step function. Since the scattering matrix is 
energy-dependent, the current I depends on the way the 
bias is applied and not just on its magnitude V . A fac- 
tor 4 has been added to the right hand side of Eq. I|15|) 
to take into account for spin and valley degeneracy. Ap- 
parently, there are terms in Eq. (|15fl that do not depend 
on time (where m = ml) and terms that oscillate as a 
function of time (where m ^= m'). 

We will now show that the dc limit of Eq. 115fl contains 
all the information of the energy dependent transmission 
of the underlying relativistic quantum dynamics. In that 
limit, only terms with m = m' survive in Eq. (|15fl and 
the differential conductance reads 
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The latter equation shows that the combination of the 
presence of the bias V and the ac signal V ac allows to 
extract the energy-dependence of the transmission coef- 
ficients (in principle) channel by channel. Typical values 
for lo that can be used experimentally in practice are up 
to 10 — 30 GHz. The bias can be controlled with any de- 
sired precision relative to ksT down to mK temperature. 

Furthermore, we find that G has an interesting feature 
which might give rise to potential applications. As illus- 
trated in Fig. |3 steps appear in dG/dV if the bias V 
equals multiples of the ac frequency uj. The steps show 
up whenever the number of propagating modes in a side 
band m vanishes. That is when the corresponding en- 
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ergy of charge carriers hits the Dirac point, where the 
density of states vanishes. The magnitude of the steps 
depends on the weight of the different side bands (given 
by the Bessel functions) and is non-universal. Never- 
theless, the appearance of these steps can be used as a 
sensitive detector for finite frequency noise, similar to the 
SIS junction in Ref. [2J}. The advantage of graphene as a 
finite frequency noise detector is that it has no frequency 
limit, whereas the SIS detector is limited by the size of 
the superconducting gap. 

At this point, we mention that our results are obtained 
using scattering theory of non-interacting Dirac fermions. 
As shown by Pedersen and Buttiker [23 (for the cor- 
responding problem based on the Schrddinger equation) 
screening of near-by gates can change the predictions of a 
scattering theory based on non-interacting particles. Al- 
though screening is expected to be reduced close to the 
Dirac point in graphene |29|. where the interesting pre- 
dictions illustrated in Fig. [21 appear, a careful analysis 
of mesoscopic ac transport in graphene in the presence 
of screening by near-by gates is an interesting research 
project by itself. 



IV. RELATION TO ZITTERBEWEGUNG 

Let us first explain the signature of ZB in the current 
in the absence of an oscillating potential (as introduced 
in Eq. (|2J)) and then show why our particular choice of 
the oscillating potential has some relevance as far as ZB 
is concerned. In the absence of the potential fi(x), the 
time evolution of the electron field operator that obeys 
Eq. can be written as 



*p(*) 



with p = (p x ,Py), p 



(17) 



*L ±} (*) 



P 2 y, 



and 



(l±(p -a)/p)% 



where a = {<7 x ,a y ) is a vector of Pauli matrices. A 
straightforward calculation of the current operator in the 
Heisenberg picture J[| shows that an oscillatory compo- 
nent in time exists due to an interference of \& p + ^ (t) and 

\Jjp \t) solutions of the DE. Therefore, in order to see 
a signature of ZB in the current, it is important that 
electron-like and hole-like solutions interfere. More gen- 
erally, an oscillatory component of the current arises if 
the electron field operator contains solutions to the DE 
at different energies. Consequently, if we calculated the 
current operator corresponding to a wave function solu- 
tion to the DE with a fixed energy e, then there would 



be no oscillatory component of the current operator left 
and, thus, no sign of ZB in the current. Importantly, this 
is the general situation in ballistic transport in graphene 
if a plane wave solution (of a particle with energy e) to 
the DE is injected from one region to another. Thus, we 
conclude that ballistic dc transport in graphene shows no 
direct signature of ZB. 

In principle, the Tien-Gordon setup of photon-assisted 
electron transport for Dirac electrons, illustrated in 
Fig. ^ can be used to generate the desired state. The 
reason is that the ac signal stimulates the absorption and 
emission of photons, which in turn lead to a population 
of different side bands around the Fermi energy. Then, 
the resulting electron field operator in region Gr (in), 
Eq. (|10l) . is precisely what is needed to observe ZB in 
graphene. However, the detection of the resulting os- 
cillating current, Eq. I|15[) . is not sufficient to prove the 
existence of ZB in graphene. The reason is that a prepa- 
ration of a state that is a superposition of different energy 
solutions to the free Schrodinger equation also yields a 
similar oscillating current (see, for instance, Eq. (16) of 
Ref. H5). The only difference to the free Schrodinger 
case as compared to the Dirac case is the peculiar energy 
dependence of the transmission in the latter. Therefore, 
the Tien-Gordon setup analyzed in this paper is poten- 
tially relevant to the detection of ZB, but our analysis in- 
dicates that the truly fundamental signature of ZB needs 
to be identified more precisely. 



V. CONCLUSIONS 

We have emphasized in this article that the energy 
dependence of the transmission (typical for mesoscopic 
graphene systems) plays a crucial role for the unexpected 
transport phenomena predicted in these devices, for re- 
cent examples see Refs. 0, 1^ E3, E3, E3, E3, E3, 
Therefore, it is desirable to directly probe the energy de- 
pendence of the transmission. We have demonstrated 
that photon-assisted transport is the natural way to do 
so. Furthermore, we have pointed out that a possible ap- 
plication of our setup is to use graphene as a spectrome- 
ter for finite frequency noise. Finally, we have discussed 
the relevance of photon-assisted transport in graphene to 
the observability of ZB. Our conclusion is that photon- 
assisted transport can be used to create the conditions 
to observe ZB. However, our analysis also shows that the 
truly fundamental signature of ZB (which is important 
for its detection) needs to be pinpointed in future studies. 
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ported by the Dutch Science Foundation NWO/FOM. 
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